Sequence co-evolutionary information is a natural partner to minimally-frustrated models of biomolecular dynamics

Experimentally derived structural constraints have been crucial to the implementation of computational models of biomolecular dynamics. For example, not only does crystallography provide essential starting points for molecular simulations but also high-resolution structures permit for parameterization of simplified models. Since the energy landscapes for proteins and other biomolecules have been shown to be minimally frustrated and therefore funneled, these structure-based models have played a major role in understanding the mechanisms governing folding and many functions of these systems. Structural information, however, may be limited in many interesting cases. Recently, the statistical analysis of residue co-evolution in families of protein sequences has provided a complementary method of discovering residue-residue contact interactions involved in functional configurations. These functional configurations are often transient and difficult to capture experimentally. Thus, co-evolutionary information can be merged with that available for experimentally characterized low free-energy structures, in order to more fully capture the true underlying biomolecular energy landscape

SMOG 2: A versatile software package for generating structure-based models

Molecular dynamics simulations with coarse-grained or simplified Hamiltonians have proven to be an effective means of capturing the functionally important long-time and large-length scale motions of proteins and RNAs. Originally developed in the context of protein folding, structure-based models (SBMs) have since been extended to probe a diverse range of biomolecular processes, spanning from protein and RNA folding to functional transitions in molecular machines. The hallmark feature of a structure-based model is that part, or all, of the potential energy function is defined by a known structure. Within this general class of models, there exist many possible variations in resolution and energetic composition. SMOG 2 is a downloadable software package that reads user-designated structural information and user-defined energy definitions, in order to produce the files necessary to use SBMs with high performance molecular dynamics packages: GROMACS and NAMD. SMOG 2 is bundled with XML-formatted template files that define commonly used SBMs, and it can process template files that are altered according to the needs of each user. This computational infrastructure also allows for experimental or bioinformatics-derived restraints or novel structural features to be included, e.g. novel ligands, prosthetic groups and post-translational/transcriptional modifications. The code and user guide can be downloaded at http://smog-server.org/smog2

Proteinlike behavior of a spin system near the transition between ferromagnet and spin glass

A simple spin system is studied as an analog for proteins. We investigate how the introduction of randomness and frustration into the system effects the designability and stability of ground state configurations. We observe that the spin system exhibits protein-like behavior in the vicinity of the transition between ferromagnet and spin glass. Our results illuminate some guiding principles in protein evolution.Comment: 12 pages, 4 figure

Large phenotype jumps in biomolecular evolution

By defining the phenotype of a biopolymer by its active three-dimensional shape, and its genotype by its primary sequence, we propose a model that predicts and characterizes the statistical distribution of a population of biopolymers with a specific phenotype, that originated from a given genotypic sequence by a single mutational event. Depending on the ratio g0 that characterizes the spread of potential energies of the mutated population with respect to temperature, three different statistical regimes have been identified. We suggest that biopolymers found in nature are in a critical regime with g0 in the range 1-6, corresponding to a broad, but not too broad, phenotypic distribution resembling a truncated Levy flight. Thus the biopolymer phenotype can be considerably modified in just a few mutations. The proposed model is in good agreement with the experimental distribution of activities determined for a population of single mutants of a group I ribozyme.Comment: to appear in Phys. Rev. E; 7 pages, 6 figures; longer discussion in VII, new fig.

Universal behavior of localization of residue fluctuations in globular proteins

Localization properties of residue fluctuations in globular proteins are studied theoretically by using the Gaussian network model. Participation ratio for each residue fluctuation mode is calculated. It is found that the relationship between participation ratio and frequency is similar for all globular proteins, indicating a universal behavior in spite of their different size, shape, and architecture.Comment: 4 pages, 3 figures. To appear in Phys. Rev.

Protein sequence and structure: Is one more fundamental than the other?

We argue that protein native state structures reside in a novel "phase" of matter which confers on proteins their many amazing characteristics. This phase arises from the common features of all globular proteins and is characterized by a sequence-independent free energy landscape with relatively few low energy minima with funnel-like character. The choice of a sequence that fits well into one of these predetermined structures facilitates rapid and cooperative folding. Our model calculations show that this novel phase facilitates the formation of an efficient route for sequence design starting from random peptides.Comment: 7 pages, 4 figures, to appear in J. Stat. Phy

Hydrophobic and ionic-interactions in bulk and confined water with implications for collapse and folding of proteins

Water and water-mediated interactions determine thermodynamic and kinetics of protein folding, protein aggregation and self-assembly in confined spaces. To obtain insights into the role of water in the context of folding problems, we describe computer simulations of a few related model systems. The dynamics of collapse of eicosane shows that upon expulsion of water the linear hydrocarbon chain adopts an ordered helical hairpin structure with 1.5 turns. The structure of dimer of eicosane molecules has two well ordered helical hairpins that are stacked perpendicular to each other. As a prelude to studying folding in confined spaces we used simulations to understand changes in hydrophobic and ionic interactions in nano droplets. Solvation of hydrophobic and charged species change drastically in nano water droplets. Hydrophobic species are localized at the boundary. The tendency of ions to be at the boundary where water density is low increases as the charge density decreases. Interaction between hydrophobic, polar, and charged residue are also profoundly altered in confined spaces. Using the results of computer simulations and accounting for loss of chain entropy upon confinement we argue and then demonstrate, using simulations in explicit water, that ordered states of generic amphiphilic peptide sequences should be stabilized in cylindrical nanopores

Molecular simulations suggest a force-dependent mechanism of vinculin activation

Focal adhesions are dynamic constructs at the leading edge of migrating cells, linking them to the extracellular matrix and enabling force sensing and transmission. The lifecycle of a focal adhesion is a highly coordinated process involving spatial and temporal variations of protein composition, interaction, and cellular tension. The assembly of focal adhesions requires the recruitment and activation of vinculin. Vinculin is present in the cytoplasm in an autoinhibited conformation in which its tail is held pincerlike by its head domains, further stabilized by two high-affinity head-tail interfaces. Vinculin has binding sites for talin and F-actin, but effective binding requires vinculin activation to release its head-tail associations. In migrating cells, it has been shown that the locations of vinculin activation coincide with areas of high cellular tension, and that the highest recorded tensions across vinculin are associated with adhesion assembly. Here, we use a structure-based model to investigate vinculin activation by talin modulated by tensile force generated by transient associations with F-actin. We show that vinculin activation may proceed from an intermediate state stabilized by partial talin-vinculin association. There is a low-force regime and a high-force regime where vinculin activation is dominated by two different pathways with distinct responses to force. Specifically, at zero or low forces, vinculin activation requires substantial destabilization of the main head-tail interface, which is rigid and undergoes very limited fluctuations, despite the other being relatively flexible. This pathway is not significantly affected by force; instead, higher forces favor an alternative pathway, which seeks to release the vinculin tail from its pincerlike head domains before destabilizing the head-tail interfaces. This pathway has a force-sensitive activation barrier and is significantly accelerated by force. Experimental data of vinculin during various stages of the focal adhesion lifecycle are consistent with the proposed force-regulated activation pathway

Atomistic simulations indicate the functional loop-to-coiled-coil transition in influenza hemagglutinin is not downhill

Influenza hemagglutinin (HA) mediates viral entry into host cells through a large-scale conformational rearrangement at low pH that leads to fusion of the viral and endosomal membranes. Crystallographic and biochemical data suggest that a loop-to-coiled-coil transition of the B-loop region of HA is important for driving this structural rearrangement. However, the microscopic picture for this proposed "spring-loaded" movement is missing. In this study, we focus on understanding the transition of the B loop and perform a set of all-atom molecular dynamics simulations of the full B-loop trimeric structure with the CHARMM36 force field. The free-energy profile constructed from our simulations describes a B loop that stably folds half of the postfusion coiled coil in tens of microseconds, but the full coiled coil is unfavorable. A buried hydrophilic residue, Thr59, is implicated in destabilizing the coiled coil. Interestingly, this conserved threonine is the only residue in the B loop that strictly differentiates between the group 1 and 2 HA molecules. Microsecond-scale constant temperature simulations revealed that kinetic traps in the structural switch of the B loop can be caused by nonnative, intramonomer, or intermonomer β-sheets. The addition of the A helix stabilized the postfusion state of the B loop, but introduced the possibility for further β-sheet structures. Overall, our results do not support a description of the B loop in group 2 HAs as a stiff spring, but, rather, it allows for more structural heterogeneity in the placement of the fusion peptides during the fusion process